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Surer et al. Reply: In the previous Comment flJJ] , Mobius 
and Richter (MR) criticize the secondary result of our Letter 
[2] because of our finite-size scaling (FSS) ansatz for the den- 
sity of states of the three-dimensional Coulomb glass (CG). 
Thus, our main result, namely that the CG does not exhibit 
a phase transition down to the lowest numerically accessible 
temperatures, remains unchallenged. This finding suggests an 
absence of a finite-temperature CG phase - an issue that has 
remained unsettled in the past 35 years and that has been in- 
dependently verified recently [3]. However, we state that the 
third result of our Letter, namely, the density of state (DOS) 
data for the random displacement (RD) version of the model, 
is incorrect. 

In contrast to the lack of conclusive results on the existence 
of a phase transition, there have been several studies attempt- 
ing to verify the estimate of Efros and Shklovskii (ES) [4] that 
the exponent of the DOS is exactly 8 = D — 1, where D is the 
space dimension. One central question of the field is whether 
numerical studies verify or contradict the ES prediction. As 
discussed in Ref. [2] there are many (often nonoverlapping) 
predictions on the exponent 8. Despite these differences, most 
publications agree that 8 is different from the ES prediction to 
some degree and that it seems to be disorder dependent. 

In an attempt to settle the controversy regarding the expo- 
nent 8, we evaluated the DOS of the CG by identifying states 
as close as possible to the exact ground state with the extremal 
optimization (EO) method [5]. Our results using EO, how- 
ever, did not deliver results qualitatively different from the 
published ones: 8 is different from the ES prediction to some 
degree and apparently disorder dependent. MR's own work 
Ja] also found these attributes for 8. 

In Ref. yfl we write that "the DOS can be fit very well 
with a form ~ _E| A " with 8 "close to the ES value" of 2 and 
not equal to 2. Our estimate differs up to 8%-10% from the 
ES prediction, whereas MR's [6] estimates are 10%-25% off. 
That MR prefer to describe this difference as a "clear devia- 
tion," whereas we used the words "close to," is an insubstan- 
tial matter of semantics. 

MR criticize our FSS form for the DOS, p(E,N) ~ 
b + a\E\ s . Here E is the energy, N = L D is the num- 
ber of sites with L the linear dimension of the sample, and 
b(N) = p(E = 0, N). To our knowledge, there have been no 
attempts to perform a clean FSS of the numerical data beyond 
the unmotivated (albeit conservative) 1/L cutoff for the pure 
power-law fitting function in Ref. 06J] . The latter corresponds 
to a step function in the FSS form, something that could be 
deemed as unphysical. This is mirrored by the large cutoff- 
dependent variations of the exponents when performing fits to 
the pure power laws [6j]. To avoid such artifacts, we chose to 
smoothly connect the low- and high-energy asymptotic behav- 
ior with the aforementioned formula. We do, however, agree 
with MR that the finite-size effects of p(E = 0, N) are in- 
deed puzzling and require further study; i.e., simulations of 


much larger systems would be needed to verify the b(N) de- 
pendence. 

Although MR's results are qualitatively similar to ours and 
bring no new insights to the problem we outline some details 
of their numerical calculations that prevent a fair comparison 
of both data sets. First, MR neglect the long-ranged nature 
of the Coulomb interactions by using the minimal distance 
method. Unphysical edge effects arise without the use of re- 
summation methods that can influence the thermodynamic be- 
havior. Second, MR propose that their data differ from ours 
because the "minimum search [we used] does not consider all 
single-particle hops." This is not correct because EO allows 
for the rearrangement of any number of particles with hops of 
arbitrary length. And indeed, we verified that our results are 
stable against single-particle rearrangements. We note that 
EO was not as efficient as we expected from studies on the 
Sherrington-Kirkpatrick spin glass [5]. 

We have also analyzed the discrepancies for the RD model. 
We reviewed our work and found an error in our code. There- 
fore, the displayed DOSs [2] correspond to disorder values 
that are ~ 10 times smaller than indicated. We have repeated 
the RD simulations with the corrected code, as well as with 
a new independent code based on single-particle hops. The 
DOSs obtained with the two methods are consistent and re- 
veal that the Wigner crystal is not as robust against disorder 
as stated in Ref. yj] . 

Summarizing, the Comment does not offer new physical 
insights. We correct our statements about the RD model in 
Ref. |2|: The Wigner crystal is not as robust to disorder as 
stated. Still, our main result of a lack of a finite-temperature 
transition in the CG remains unchallenged. 
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